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Abstract. We survey known results and present estimates and conjectures for the next- 
order term in the asymptotics of the optimal logarithmic energy and Riesz s-energy of N 
points on the unit sphere in E'^+^, d>l. The conjectures are based on analytic continu- 
ation assumptions (with respect to s) for the coefficients in the asymptotic expansion (as 
N ^ oo) of the optimal s-energy. 



1. Introduction 

Let S*^ denote the unit sphere in the Euchdean space M''"'"^, d > 1. The discrete logarith- 
mic energy problem on S*^ is concerned with investigating the properties of A^-point systems 
x^, . . . , on S'^ maximizing the product of all mutual pairwise Euclidean distances 

M(xi, . . . ,XAr):=]^ |Xj - Xfcl = Y\_ \^j-^kf^ 
j=/=k l<j<k<N 

or equivalently, minimizing 

(1) Eiog(xi, . . . ,XAr):= - logM(xi, . . . ,XAr) = ^log 

over all A^-point configurations (xi, . . . , xat) on S*^. 

One goal of this article is to investigate the asymptotic expansion as goes to infinity 
of the minimum (N -point) logarithmic energy of S'^ 

fiog(§^; N):=mm {Eiog(xi, . . . , x^v) I Xi, . . . , x^ G S'^} = Eiog(x*, . . . , x^), N >2. 

We remark that in his list of "Mathematical problems for the next century" Smale [121 
Us] posed as Problem #7 the challenge to design a fast (polynomial time) algorithm for 
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generating "nearly" optimal logarithmic energy points on the unit sphere in M.^ that satisfy 
(2) Eiog(xi, . . . ,XAr) — £^iog(§^; A^) < clog for some universal constant c. 

This problem emerged from computational complexity theory (cf. Shub and Smale [H] 



The right-hand side of ([T]) is referred to as the discrete logarithmic energy of the normal- 
ized counting measure /i[xi, . . . , x^r], which places the charge at each point xi, . . . , x^r. 
By the continuous logarithmic energy of a (Borel) probability measure n supported on S'^ 
we mean 

Iiog[fA-= Jj log |x-y| d/x(x)d/i(y). 

Note that Xiog[/i] = +00 for any discrete measure fi. Classical potential theory yields that 
Xiog[/i] is uniquely minimized by the normalized surface area measure a^, Jg^ dcr^ = 1, over 
the class of all (Borel) probability measures /i supported on S'^; that is, 

(3) V,,,{S'y.= inf {Tiog[/.] I /i G A^(§^)} = Xiog[a,] = log ^ + ^ [m - ^{d/2)] , 

where ^{z):=T' (z) / T(z) is the digamma function. The last expression in (|3| follows, for 
example, from Eq. (2.26)]. 

It is known that the minimum A^-point logarithmic energy of S'^ satisfies 

(4) V rs^l 1 logiV ci giog(S'^;iV) , 1 log AT c, 

for some constants ci, C2 depending on d only. The lower bound follows from |39]. The 
upper bound follows from an averaging argument as in the proof of Theorem 1 in j27] . 
which is based on equal area partitions [SU]. These bounds give the correct form of the 
second-order term in the asymptotics of the minimum A^-point logarithmic energy of S^. It 
should be mentioned that |3B] also cites [221 P- 150] as a reference to a more general method 
to obtain lower bounds for minimum logarithmic energy valid over more general Riemann 
surfaces, which in the case of the 2-sphere also gives the lower bound in ^ (cf. [Ml Sec. 3]). 
In a recent paper the first author improved the lower bound for higher- dimensional spheres. 
In [T^ Lemma 1.1] it was shown that 

(5) 5^>^„,(svi!^-f + a(iv-n n^o., 

where the positive constant does not depend on A^ and is given by 

C'-V I md)T{l+[d/2\-d/2) l^^l^. 



the number e satisfies < e < 1 {d even) or0<£<l/2((i odd). Combining the upper 
bound in ^ and the lower bound in ^ we obtain the following asymptotic expansion of 
the minimum A^-point logarithmic energy of S'': 

(6) £iosi^'';N) = ViogiS'')N^ -^N\ogN + OiN), N ^ 00. 



THE NEXT-ORDER TERM 3 

We remark that the minimum A^-point logarithmic energy of the unit circle § is attained 
at the A^-th roots of unity and (cf., for example, |15j ) 

Siog{S;N) = -NhgN, N>2. 

(In this case Viog(S) = 0.) 
In view of ^ it is tempting to ask if the following limit exists: 



(7) lim 4 



£iog(§''; N) - yiog(§'') + log N 



In particular, it is plausible and consistent with the lower bound ([s]) that this limit is 
negative for all c? > 2 if it exists. Indeed, in [HSl IHU] it is shown (also see below) that on the 
unit sphere in M.^ (that is, d = 2) the square-bracketed expression in ([T]) can be bounded 
by negative quantities from below and above. 

In formulating our Conjecture |4] for £^iog(§^; A^) (see Section |4]), we will make heavy use 
of the observation that the discrete logarithmic energy is the limiting case (as s — > 0) of 
the Riesz s- energy 

E,(xi,...,X7v):= ^ = 2 ^, 

by means of (1/|-|'^ — — i- log(l/|-|) as s — )■ 0. For A^ > 2, we shall also consider the 
optimal N -point Riesz s- energy of a compact set A C M'^^^ defined by 



(8) Ss{A;Ny.-- 



min {Es(xi, . . . , x^r) | Xi, . . . , xjv ^ A} for s > or s = log, 
max {Es(xi, . . . , xjy) | Xi, . . . , xjv G ^} for s < 0, 



where we note that £o{A; N) = N"^ — N, which is attained by any A^-point configuration 
on A. Furthermore, it is known (cf. [12]) that 



(9) ^SsiAN) 

as 



= £iog(A,Ar) (Ar>2). 

s=0+ 

As in the logarithmic case, the optimal continuous s-energy of S'^ for -2< s <d IS given 

by 

(10) V.(S^):^ //^d..(x)dMy)^2--r((''+^'/^'^((<'-"/^' 



|x-y|" ^/^TT{d-s/2) 

and the s-capacity of S'^ is given by cap^(S'^):=l/ys(§'^). 

In Sections [3] and |4] we will discuss estimates and asymptotics for the optimal A^-point 
s-energy on S'^. 

2. The logarithmic energy on the unit sphere in 
Let Aiog(A^) A^ be the remainder term in 

£iog(§'; N) = \4og(S') - ^N\ogN + Aiog(Ar) N. 
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Rakhmanov, Saff and Zhou [3Sl EH] proved the following estimate^ 

-0.22553754... 



(11) 
(12) 
where 



liminfA,og(iV)>-ilog (l-e"^)"' 



TT 



1 1 7r-\/3 
limsup Aiog(iV) < - - log — - + — = 
7V-s>oo ^ ^ 4vo 
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and they stated the following conjecture: 



Conjecture 1 ([38j). There exist constants Ciog2 (^iT'd D\og2, independent of N, such that 
(13) 

£iog(§'; iV) = Q - log 2^ - iiVlogiV + Ciog,2iV + Aog,2 logiV + 0(1) as iV oo. 

The authors of |38j supported their conjecture by fitting the conjectured formula to 
the data obtained by high-precision computer experiments to find extremal configurations 
(A^ < 200) and their logarithmic energy by minimizing the absolute £i-deviation. This 
leads to the following approximatior|^ 



^iog(S';iV) 



^-log2j N^-^NlogN 



0.052844iV + 0.27644. 



(In [38] it is remarked that the logarithmic term was ignored in the fitting algorithm 
because it did not appear to be significant.) 

So far no conjecture for the precise value of Ciog,2 in (13) has been given. Conjecture |4 



in Section |4] fills this gap. Our approach for arriving at this conjecture is as follows: Since 
the logarithmic energy of an iV-point configuration X^r can be obtained by taking the 
derivative with respect to s of the Riesz s-energy of and letting s go to zero, we shall 
differentiate the asymptotic expansion of the minimal N -point Riesz s-energy of §^ (see 



and then let s go to zero from the right to derive plausible asymptotic formulas for 
the X-point logarithmic energy of S^. In this derivation there is an implicit interchange of 
differentiation and minimization which is formally justified by (|9]). 

In order to illustrate this approach, we begin with the simplest case, which is the unit 
circle §. Recently, we obtained in [IS] the complete asymptotic expansion of the Riesz 
s-energy Cs{N) of N-th roots of unity, which are the only (up to rotation) optimal energy 
points on the unit circle if s > —2, s 7^ 0, as well as for the logarithmic case. In fact, the 



*We remark that there is a typographical error in the sign of the second ratio in the upper bound in 
[3H1 Eq. (3.10)]. 

''"Note that the logarithmic energy in [551 EH] is precisely half of ours. 
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expansion is valid for s G C with s not zero or an odd positive integer: 

(14) (27r)^ {27iy ^ 

+ a,p(iV"'+'^''~'^) as AT ^ cx) (p = 1, 2, 3, . . . ), 

where ({s) is the classical Riemann zeta function. The coefficients an{s), n > 0, satisfy 
the generating function relation 

{-irBiiis/2) 



smvr^ 



nz 



E 

n=0 



a„{s z 



2n 



\z\ < 1, s G C, 



a„[s 



{2n) 



i2n 



,2n 



Here, B^\x) denotes the so-called generalized Bernoulli polynomial of degree n. In par- 
ticular, one has the following recurrence relation for n > 1 (cf. Luke [31, p. 34f]) 

n-l 



(15) 



m=0 



2n-l \ B. 



^2m+2 jj{2p) 



2m + 11 2m + 2 



2n-2-2m 



where Bq = 1, Bi = —1/2, B2 = 1/Q, B^ = 0, . . . are the Bernoulli numbers. The 
constant Vs{S) is the Riesz s-energy of the unit circle if — 2 < s < 1, s 7^ 0, and its analytic 
continuation to C \ {1, 3, 5, 7, ... } otherwise: 

2-r((i-s)/2) 



(16) 



VsiS) 



-, s G C, s 7^ 1, 3, 5, 



1. 



/^r(i - s/2) 

Using properties of the digamma function and the Riemann zeta function, we obtain 

'1 



ds 



s-s>0 



2C(^) 
ds { (27r)^ 



l+s 



{27ry 



{^{l-s/2)-^{{l-s)/2))-\og2 
Cis 



2 



logN 



Cis) 



log(27r) 



= 0, 

s-s>0 

-AT log N. 



Moreover, using that the Riemann zeta function vanishes at negative even integers and 



5^J)(0) = for A; > 1 (cf. (pi)), we obtai 



cam 



2a J s) 



N 



l+s-2n 



Cis - 2n) 



ds [ (27r 

Thus, putting everything together, we arrive at 

dCsiN) 



for n = 1, 2, 3, ... . 



^iog(S; N) 



ds 



-N\ogN, 



as we expected. 

A second fundamental observation is that the coefficients of the powers of A^ in the 
asymptotic expansion of CsiN) can be extended to meromorphic functions on C. In this 
way one sees that the leading two terms "swap places" when moving s from the interval 
(—2, 1) to interval (1, 3). (Note that we assume that CgiN) = N{N — 1) at s = 0.) Thus, 
each term's singularity at s = 1 can be avoided by moving away from the real line. In other 
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words, the second term of the asymptotics of Cs{N) in the finite energy case s G (—2, 1) is 
the analytic continuation of the leading term in the hyper-singular case s G (1, 3) and vice 
versa. 

This is the motivation for trying the same ideas for the rf-sphere in Section |4j 

3. Asymptotics of optimal Riesz s-energy 

3.1. The dominant term. The leading term of the asymptotic expansion (as — )■ oo) 
of the maximal (if s < 0) and minimal (if s > 0) Riesz s-energy of a compact set A in MP is 
well understood if A has finite s-energy (that is, positive s-capacity). For sets of vanishing 
s-capacity the leading term is rather well understood in the sense that the existence of the 
coefficient has been established for a large class of compact sets by the second and third 
authors, but the determination of this coefficient ioT s > d (except for one- dimensional 
sets) has to date remained a challenging open problem. 

The general theory for (the energy integral associated with) the continuous s-potential 
{s < 0) covered in Bjorck [9] provides Frostman-type results, existence and uniqueness 
results for the equilibrium measure Hs, and characterization of the support of fis for general 
compact sets A in W. Of particular interest is the observation that the support of fis is 
concentrated in the extreme points of the convex hull of A if s < —1 and for s < —2 any 
maximal distribution (there is no unique equilibrium measure anymore for s < — 2) consists 
of no more than {p + 1) point masses. The singular Riesz s-potential on S'^, < s < ci, is 
the subject of classical potential theory (see, for example, Landkof [28j) with the value of 
the energy integral for S*^ being already studied by Polya and Szego [33] . The range of the 
potential theoretical regime is thus limited to —2 < s < d. 

The first results for the (i-sphere in the hypersingular case s > d were given in [27]. This 
included the leading term for the exceptional case s = d, bounds for the leading term for 
s > d, separation results, and conjectures for the leading coefficient for the 2-sphere. In 
subsequent work the existence of the leading term for s > d was proven for the class of 
(i-rectifiable sets ([2ll|25]) as well as for weighted Riesz s-energy on such sets (pT]). 

The potential-theoretic regime —2 < s < d. A standard argument from classical potential 
theory involving the concept of transfinite diameter yields that the positive N-th gener- 
alized diameter o/ S"^, Dj^{E)'^):=N{N — l)/£^s(S'^; A^), forms a monotonically decreasing 
(increasing) sequenc^ bounded from below (above) if s > (s < 0). This implies the 
existence of the limit D*:= lim7v-!.oo -D|^r(S'^) called the generalized transfinite diameter of 
S'^ introduced by Polya and Szego in ^35J, where it is further shown that the s-capacity 
and generalized transfinite diameter are equal. 

Theorem 1. Let d > 1. Then for —2<s<d 

lim Ss{S'^]N)/N^ = VsiS"^) = l/cap,(§'^) = 1/D'. 



*This monotonicity holds more generally for any compact set A. 
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Thus, in the potential-theoretic regime the dominant term grows hke A^^ as — ?■ oo 
and its coefficient in the asymptotic expansion of the optimal s-energy of S'^ encodes the 
s-capacity as well as the generalized transfinite diameter of S'^. 

For future reference, we remark that V"s(§'^) can be analytically extended to the complex 
s-plane except at the simple poles at s = d + 2k, A; = 0, 1, 2, ... if d is odd and for 
/c = 0, . . . , I — 1 if (i is even. The residue at s = d + 2k, regardless of the parity of d, is 

(17) (_i)fc2-2'=-i r((rf + l)/2)/ [V^k\ Tid/2 - k)] . 

We will denote this meromorphic extension with the same symbol Vs(§'^). 

The hypersingular case s > d. Since foT s > d the s-energy integral for every positive Borel 
probability measure supported on S"^ is +oo, potential theory fails to work. (The boundary 
or exceptional case s = d can still be treated using a particular normalization of the energy 
integral and a limit process as s approaches d from below, see [16]). 

The dominant term of the asymptotic expansion of the minimal s-energy grows like 
A^^ log in the boundary case s = d. The coefficient is also known. 

Theorem 2 (boundary case s = d, [2Z])- 

...^ qs-^; AT) _ ?/,(B'^) _ i r((d+l)/2) 

^ ' N^oo m\ogN Hdi^'^) d ^V{d/2) ' 

Here and hereafter, T-Ldi') denotes ci-dimensional Hausdorff measure in W, p > d, normal- 
ized such that a c?-sided cube with side length 1 has "Hrf-measure equal to 1. In particular, 
'Hrf(B'^) denotes the volume of the unit ball in M*^ while 'Hd{^'^) denotes the surface area of 
the unit sphere in M*^"*"^. 

In [27] the order of the growth rate of Ss{Ei'^; N) was established: there are constants 
Ci,C2 > such that 

(19) Ci A^i+^Z'^ < £,(§'^; N) < C2 A^i+^Z'^ {s>d>2). 

The second and third author [25] showed that the limit of the sequence Ss{S^; N)/N^~^'^^'^, 
indeed, exists. More generally, the following result holds, which has been referred to as the 
Poppy-seed Bagel Theorem because of its interpretation for distributing points on a torus: 

Theorem 3 ([23 [II])- Let d > 1 and A G W an infinite compact d-rectifiable set. Then 
for s > d 

(20) hm Ss{A; N)/N'+^/' = Cs,d/ [HdiA)]'/' , 
where Cs^d is a finite positive constant (independent of A). 

By d-rectifiable set we mean the Lipschitz image of a bounded set in M"^. 
In particular. 



s > d. 
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In Thm. 3.1] it is shown that Cg^i = 2(^(s). For d > 2 the precise value of Cs^d is 
not known. The significance (and difficulty of determining Cg^d) is deeply rooted in the 
connection to densest sphere packings. Let 6n denote the best-packing distance of A^-point 
configuration on S'^. It is shown in pilj that 

(21) hm [C,,,]i/^ = l/Coo,d, C^,d-= lim iW'^^^ = 2 [A,/Hd(B'^)]i/^ 

s— >oo N—^co 

where is the largest sphere packing density in M'', which is only known for three cases: 
Ai = 1, A2 = 7i/^/l2 (Thue in 1892 and L. Fejes Toth [2T]), and A3 = n/^/l8 (Kepler 
conjecture proved by Hales ^23j). 

It is not difficult to see that the Epstein zeta function Ca{s) for a lattice A C M'^ defined 
for s > d by 

Ca(.) := J2 

yields an upper bound for Cs^ as we now show. Let Q denote a fundamental parallelotope 
for A. For n G N, the intersection of Q and the scaled lattice (l/n)A contains exactly 
N = n'^ points and, for s > d, we have 

E,ix^)=j2 E E \^-yr = n'^XAis) = N'^-^/''CA{s). 

xeX„x^yeX„ xeX„ x^y6{l/n)A 



Referring to (20) with A = Q, one obtains the following result. 
Proposition 1. For s > d, 

(22) C,,d<min|Ar/'^CA(s), 

A 

where the minimum is taken over all lattices A C M*^ with covolume |A| > 0. 
In particular, as shown in 




r ^ \ ^ - \ \ f , A r gs(§^iV) ^ (v/3/2)-/^Ca.(^) 
Cs,2 < I — 1 Ca.(^) and Jim ^^^^ < , . > 2, 

where A2 is the hexagonal lattice consisting of points of the form m(0, 1) + n(l/2, V3/2) 
with m and n integers. 

For most values of d we do not expect equality to hold in (22); that is, we do not expect 
lattice packings to be optimal (especially for d large where it is expected that best packings 
are highly 'disordered' and far from being lattice arrangements, cf. |47]); however, recent 
results and conjectures of Cohn and Elkies [H], Cohn and Kumar [19] and Cohn, Kumar, 



and Schiirmann pO] suggest that equality holds in (22) for d = 2,4,8 and 24 leading to 
the following: 

Conjecture 2. For d = 2, 4, 8 and 24, Cg^d = \-^d\^^'^CAd{s) for s > d, where denotes, 
respectively, the hexagonal lattice, D4, Eg, and the Leech lattice. 

For d = 2, this conjecture appears in 
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3.2. The second-order term. The only known resuhs so far are estimates of the differ- 
ence £^s(S^; A^) — Vs{S'^)N'^ in the potential-theoretic regime. At the end of this section we 
present a lower bound for the optimal d-energj, the first hypersingular case. 

The average distance problem (s = —1) on a sphere was studied by Alexander Oil] (d = 
2) and Stolarsky [l5] {d > 2, see citations therein for earlier work), later by Harman [26] 
{d > 2, correct order and signs of the bounds up to logA^ factor in upper bound) and 
Beck [7J {d > 2, settled the correct order and signs of the bounds); generalized sums of 
distances (—2 < s < 0) were studied by Stolarsky |44j on rf-spheres and by Alexander and 
Stolarsky [3] , Alexander [3] for general compact sets in W and general energy functionals; 
eventually, Wagner ([IS] upper bound and [IH] lower bound) arrived at 

-a A^^+"/°' < SsiS'^; N) - VsiS'^) < -(3 N^^'l'^ (-2 < s < 0). 

where a and (3 are positive constants depending on s and d but not A^. 

For the singular Riesz s-potential on Wagner [19] found the upper bound 

A^) < K(S^) A^^ - C2 N^+'/^, > (0 < s < 2). 

The method of the alternative proof in 1^38] was generalized in [27] leading to 
^(S'^; A^) < VsiS'^) - C'2 N^+'/^, > {0<s<d). 

Wagner [IH] also gave the lower bounds 

'CiA^i+^Z'^ d-2<s<d, 



£si^ ■,N) > VsiS'')N' - 



C[ A^i+^/(^+2) 0<s<d-2, d>3, 



which were improved by the first author [13j. All results combined together lead to the 
correct order of growth for the second-order term: 

Proposition 2. Let d > 2. Then for each —2<s<d, there exist constants c, C > 
which depend on s and d, such that 

-cN'^+'l'^ < SsiS'^; N) - K(§'^) A^^ < -CN^+'Z"^, N>2. 

Next, we present bounds for the hypersingular case s > d which follow from a careful 
inspection of the proof of the dominant term in [27] . See also [13] for the potential-theoretic 
case < s < d. 

Proposition 3. Let d > 2. Then, as N ^ 00, 

-c{d)N'+OiN'-'/'\ogN) < £di^'-^N)-^^NHogN < ^^^^ NHog\ogN+0{N') 
where the constant c{d) is given by 

<*-|g|{l->o.|g.<*l*W2)-W-.og2]}>0. 

(Recall that ip denotes the digamma function.) 
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Remark. For d = 2 one has c(2) = 1/4 and an 0{N) term in the lower bound instead of 
0{N log N). 

The proof of Proposition |3] is given in Section [6] along with proofs of other new results 
stated in this section. 

Following the approach leading to Proposition [3] we obtain for s > ci the following crude 
estimate, which, curiously, reproduces the conjectured second term but only provides a 
lower bound for the leading term (that is, for the constant Cg^d in the leading term). 

Proposition 4. Let d>2 and s > d. Then, for {s — d)/2 not an integer, 

SsiS"^; N) > As,d N^^'''^ + Vs{^'^) + 0{N^+''^^^''^) as N oo, 

where 



A.. 



d 



s — d 



ir((rf + i)/2)r(i + (s-rf)/2) 

2 v^r(l + s/2) 



s/d 



Note that for s > + 2 in the above proposition the 0{N^^^^^'^^^''-) term dominates the 
A^^-term. 

The upper bound of ^^s(S'^; N) in the hypersingular case s > c? is in the spirit of Propo- 
sition m 

Proposition 5. Let d>2 and s > d. Then, for {s — d)/2 not an integer. 



s/d 



['Hd{S^){l-d/s)_ 

The constant Ag^d in Proposition |4] and the corresponding constant in Proposition [5] 
provide lower and upper bounds for the constant Cs^d appearing in Proposition [2j 

Corollary 4. Let d>2 and s > d. Then for {s — d)/2 not an integer 



A,,d < 



< 



Hd{S'i){l-d/s) 



s/d 



[Hd{S<^)Y^' 
where Ag^^ is given in Proposition^ 

For d = 2, the above bounds for Cs,d/[Hdi^^)Y^'^ reduce to 

(23) 2-^/^5-^/7 [^/2 - 1] < Cs,2/ [47r]^/' < 2'' [s/ {s - 2)]^/' 

4. Conjectures for the Riesz s-energy 



A straightforward generalization of the asymptotics for the unit circle ( 14 ) would be 

p 



(24) 



C. 



s,d 



[Udmr" 

+ 0,^d,p{N^^"''^~^^''^~'^''^) 



^l+s/d ^ ^^^^^ ^■^^i+s/d-2n/d 

n=l 

as N 00 
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for s not a pole of KlS'^) given in (10) and Cg^d as defined in (20). The exceptional cases 
are caused by the simple poles of as a complex function in s. As s — d, one of the 

terms in the finite sum has to compensate for the pole of the A^^ term, thus introducing a 
logarithmic term (cf. the motivation for Conjecture [s] in Section [?]). Recall that for even 
d there are only finitely many poles suggesting that there might only be finitely many 
exceptional cases for even d. 

It is possible that the higher-order terms (p > 1) may not exist in this formj^ the 
presence of the A^^ and A^^+^Z*^ is motivated by known results and the principle of analytic 
continuation. In the potential theoretic case —2 < s < d, the leading term is the A^^- 
term and KlS'^) equals the continuous s-energy of S*^ (see Theorem [l| and is its analytic 
continuation elsewhere. By the same token, the leading term in the hypersingular case 
s > d is the A^^+'^/'^-term (see Theorem 2| whose coefficient Cs^d is assumed to have an 
analytic continuation to the complex s-plane as well. 

Conjecture 3. Let d > 2. For —2 < s < d + 2, s ^ d, there is a constant Cg^ such that 

CsA 



(25) 8s{^'^;N) = Vs{^'^)N^ + 



N^+^/d + o(Ari+"/'^) as N ^oo, 



where, for s > d, the constant Cs^ is the same as that in (20). Furthermore, ford = 2, 4, 8 
and 24, Cg^ = \Ad\^^^CAj^{s) , where Ad is as in Conjecture For d = 2, — 2 < s < 4, and 
s ^2, (25) reduces to 



(26) £.(§'; AT) 



3/2r^'cA,(^) 

(4vr: 



s/2 



Ari+-/2 + o{N^+''^) asN-^oo. 



Remark. Note that s = in (25) does not refer to the logarithmic case but rather to 
£o(S'^; N) = N{N - 1), from which we deduce that Co,, 



-1. 



For s in the range —2 <s<d + 2(s^d), the A^^-term and N^'^'^^'^-teim are the two 
leading terms, interchanging their role as dominant term as s passes by the boundary case 
s = d, where V^(S'^) has a pole. Its next pole ioi s > d occurs at s = ci + 2. For higher 
values of s {s > d + 2) the A^^-term is conjectured to be dominated by other powers of A^. 

Remark. Starting with the assumption that the Riesz s-energy of A^ points is approximately 
given by A^ times the potential $ (or "point energy" ) created by all other A^ — 1 points at 
a given point and using a semicontinuum approximatioiij^] to approximate Berezin [8] 



^In fact, Michael Kiessling from Rutgers University pointed out that higher-order terms might not exist, 
instead one could have an oscillating term as N grows. Moreover, numerical results in Melnyk et al [55] 
suggest that for N fixed there might be an abrupt change of optimal configuration when an increasing 
s passes certain critical values that depend on N. This behavior might also influence the existence of 
higher-order terms. 

^This technique can be found in old papers addressing problems in solid state physics, cf. p. 188f]. 
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arrived at the plausible asymptotics 

2i-s 



2-s 



1 - {n/Nf-'^^ 



X 



6 

1^ 



6 



+ N 



12 



s/2 



6 



Stt 



6 6 

^ 3^ ^ {2^y 



+ 



12 



13)^ 



+ 



The denominators [without the power] are the first 7 distances in the hexagonal lattice A2 
and the numerators give the number of nearest neighbors with the corresponding distance. 
Thus, $ is approximated by the flat hexagonal lattice using the nearest neighbors up to 
level 7 and the remaining N — n points are approximated by a continuum), where n is 
one plus the number of at most 7-th nearest neighbors. The • • • indicates that the local 
approximation can be extended to include more nearest neighbors. In fact, for s > 2 
the square-bracketed expression is the truncated zeta function for the hexagonal lattice. 
Figure [T] shows the distance distribution function for numerical approximation of a local 
optimal 900-point 1-energy configuration (cf. Womersley [51]). The first few rescaled 
lattice distances are superimposed over this graph. Note the remarkable coincidences with 
the peaks of the distance distribution function. 

The zeta function Ca2('^) appears in number theory as the zeta function of the imaginary 
quadratic field Q(^/— 3), whose integers can be identified with the hexagonal lattice A2. It 
is known (cf., for example, |T7l Ch. X, Sec. 7]) that CAal-^) admits a factorization 

(27) Ca,(s) = 6C(s/2)L_3(s/2), Res > 2, 

into a product of the Riemann zeta function C, and the first negative primitive Dirichlet 
L-Series 



(28) 



L_3(s):=l 



1111 

2s ~^ 4s 5s 7s 



Res > 1. 



The Dirichlet L-series above can be also expressed in terms of the Hurwitz zeta function 
C(s,a) at rational values a, cf. ^27j. That is. 



(29) 



L_3(5) = 3-^ [C(s, 1/3) - 2/3)] , C(^, «):= 



k=0 



{k + a) 



Res > 1 



Remark. It is understood that in (26) is the meromorphic extension to C of the right- 
hand side of (27). Since ({s) is negative on the interval [—1, 1), has a pole at s = 1, and is 
positive on (1, 00) and the Dirichlet L-Series is positive on the interval (—1, ooj^ (cf., Eq. 
(27)), it follows that Cs^2 in (26) would be negative for —2 < s < 2 and positive for s > 2. 



Based on the motivating discussion in Section 7, we propose the following conjecture for 
the logarithmic energy. 



llNote that L_3(l - 2m) = for m = 1, 2, 3, 4, . . . , cf. [SD]. 
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Figure 1. The histogram (counting number of occurrences) and the em- 
piric distribution function of the first 100000 distances of a local optimal 
900-point 1-energy configuration. The vertical lines denote the hexagonal 
lattice distances adjusted such that the smallest distance coincidences with 
the best-packing distance. 



Conjecture 4. For d = 2,A, 8, and 24, 



(30) ^iog(§^; N) = Mog(§'') N^--NlogN + Qog.diV + o{N), as N ^ oo, 



where the constant of the N-term is given by 

1 



(31) 



aog,. = ^log(H.(§'^)/|A.|) + C,(0). 



For the case d = 2, (31) reduces to 



Q„,, = 21og2+ilog^+31og^^^^3j 



-0,05560530494339251850 , . . < 0. 
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Remark. We expect, more generally, that (30) holds for arbitrary d provided that Cg^d is 
differentiable at s = 0, in which case (31) becomes 

d 



log,ii 



ds 



C. 



+ {l/d) log?/,(§'^). 



s=0 



Regarding the bounds (11) and (12), note that Ciog,2 is closer to the upper bound given 
in (12), which gives rise to the question if the related argument can be improved to give 
the precise value. 

4.1. The boundary case s = d. As with the unit circle, we expect to obtain the asymp- 
totics of the optimal Riesz energy in the singular case from the corresponding asymptotics 
of the Riesz s-energy, s ^ d and s sufficiently close to by means of a limit process s ^ d. 
This approach leads to the next conjectures whose motivating analysis is given in Section 
7. 



Conjecture 5. Let d>l. Then 



-Ham 



log N + Cd,dN' + 0{l) asN^oo, 



where 



C, 



For d = 2, this becomes 
1 



C, 



2,2 



4 



7 



- log(2v^7r) 



477 



lim 



s,d 



[Hdm] 



[7i(2/3)-7i(l/3)] 



-0.08576841030090248365 • • ■ < 0. 



Here, 7 is the Euler-Mascheroni constant and '^n{oi) is the generalized Stieltjes constant 
appearing as the coefficient jnio^/nl of (1 — s)" in the Laurent series expansion of the 
Hurwitz zeta function ({s, a) about s = 1. 

5. Numerical Results 

Rob Womersley from UNSW kindly provided numerical data, which we used to test our 
conjectures. For the logarithmic and the Coulomb cases (s = 1) on §^ results for small 
numbers {N = 4, . . . , 500) and for large numbers of points {N = (n + 1)^ points, N up to 
22801) are given. The reader is cautioned that these numerical data represent approximate 
optimal energies which we denote by £^iog(S^; A^) or ^s(§^; N). 

A general observation is the slow convergence of the sequence of s-energy values. 

5.1. Logarithmic case. Figures |2] and [3] show the convergence to the conjectured coeffi- 
cient of the iV-term (see Conjecture |4]) 



(32) 



log I 



?2. 



N) 



1 



NlogN 



/N. 



The horizontal line indicates the value of Ciog,2 given in Conjecture |4| 
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0.00 - 



-0.01 



-0.02 



-0.03 - 



-0.04 - 



-0.05 




Figure 2. The logarithmic case on = 4, 5, 6..., 500 points. The 



horizontal axis shows iV on a logarithmic scale and the quantity (32) on the 
vertical axis. The horizontal line shows the conjectured limit C\og2- 



-0.01 



-0.02 - 



-0.03 - 



-0.04 - 



-0.05 - 




50 100 



500 1000 



5000 ixlO'* 



Figure 3. The logarithmic case on §^ with iV = 4, 9, 16 . . . , 10201 points. 
The same quantities are shown as in Figure |2] 



5.2. The Coulomb case s = 1. Figures |4] and 
coefficient of the A^^"^^/^-term (see Conjecture^) 



5] show the convergence to the conjectured 



(33) 
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The horizontal hne indicates the value of Ci^2/V^- 



-1.080 



-1.085 



-1.090 



-1.095 



-1.100 



-1.105 




10 



20 



50 



100 



200 



500 



Figure 4. The Coulomb case (s = 1) on = 4, 5, 6 . . . , 500 points. 

The horizontal axis shows on a logarithmic scale and the quantity (33) on 



the vertical axis. The horizontal line shows the conjectured limit Ci^2/v4vr. 



-1.085 



-1.090 



-1.095 



-1.100 



-1.105 




50 100 



500 1000 5000 1 X 10'' 

2 



Figure 5. The Coulomb case (s = 1) on §^ = 4, 9, 16 . . . , 22801 points. 
The same quantities are shown as in Figure |4j 
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5.3. The boundary case s = d = 2. Figure [6] shows the convergence to the conjectured 
coefficient of the A^^-term (see Conjecture ^ 



(34) |^2(S'; N)-^N^ log ivj /N^ 

The horizontal line indicates the value of 6*22- 




Figure 6. The boundary case (s = 2) on S^, iV = 4, 5, 6 . . . , 500 points. 
The horizontal axis shows on a linear scale and the quantity (34) on the 
vertical axis. The horizontal line shows the conjectured limit 6*2,2- 



6. Proofs 

In the following we set 

^^^^ ) = r((rf + i)/2)- 

Proof of lower hound in Proposition This proof follows the first part of the proof of The- 
orem 3 in • By an idea of Wagner, the (hyper) singular Riesz d-kernel is approxi- 
mated by the smaller continuous kernel l/{e + r'^Y^'^ (e > 0). Then, for x, y G S'', we have 
l/(£ + |x-y|2)'^/2 = ir^((x,y)) where Ke{t):={2-2t + e)-'^''^ . Note that Ke is positive defi- 
nite in the sense of Schoenberg [ID]; that is, it has an expansion K^if) = J2'^=o (^n{£) Pni't) in 
terms of ultraspherical (normalized Gegenbauer polynomials -P^(t) = C^(t)/C^(l), where 
X = {d — l)/2) with positive coefficients an{e) {n > 1) giving rise to the estimates 

(35) jjLk [e + |xj - Xfcl J j-fc j 

>ao{e) N^-K,{1)N, 
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where we used the positivity of i^^ and -Pq (t) = 1 (cf. Sec. 3 of [27]). Using the integral 
representation [H Eq. 15.6.1] of the regularized Gauss hypergeometric function we obtain 



ao(e 



(4 + e 



,-d/2 ryd-l^d-l 



^d Jo 



d/2-1 



l-U 



t 

d/2-1 



-d/2 



du 



= (4 + e)-'" 2'^-^'^V{d/2f ( ^^2. _ 
Ud \ a 1 + e/AJ 

Assuming that e < 2, we apply the linear transformation [H Eq. 15.8.11] with the under- 
standing that V((i/2 - k)/T{d/2 - k) is interpreted as {-lf-'^/'^+^{k - d/2)\ if d/2 - k is 
a non-positive integer: 



ao{e) 



'-T{d/2)J2 



id/2), 



2 ' f^^k\k\T{d/2-k) KA 



log - - 2 ^{k + 1) + ^{d/2 + k) + ^{d/2 - k) 



2 ujd 



^(-log^ .F, (1-^/2,^/2. 

oo 

- r(rf/2) ^ 



2 



fc=0 



k\k\T{d/2 - k) \A 



[i){d/2 + k)+ i){d/2 -k)-2^{k + 1)] . 



Note that the (non-regularized) Gauss hypergeometric function is a polynomial of degree 
d/2 — 1 if (i is even and reduces to 1 if li = 2. Using the series representation [1^ Eq. 15.2.1] 
we arrive at 

ao(e) = \—{-\oge) + — log2 + 0{e) = \ \og{A/e) + 0{e) as e ^ if d = 2, 

2 UJ2 U}2 4 



-(-log£)-^^[^(d/2)-^(l)-log2] + C(£log(l/£)) ase^0ifd>3. 

2 uJd ^d 



Using the substitution e = a^N '^1'^ and noting that -^'^(l) = s '^Z^, there follows from (35) 
(one has C(A^2-2/d^ as ^ cx) if = 2) 



1 log - /(2; a) + 0{N), 



d = 2, 

^ N^logN - f{d]a)N^ + 0{N^-^/'^\ogN), d>3, 
as — )■ oo, where for each positive integer d the function 



fid; a):-- 



[loga + ^{d/2) - ^(1) - log 2] + a- 



a > 0, 



has a single global minimum at a* = [cud-i/ idud)] ^^'^ in the interval (0, oo) with value 



fid; a*) 



1 (^d- 



d UJd 



1 - log(^^) + d{^{d/2) - ^(1) - log 2) 
d Ud 
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Let Fci:=uJd~i/ {d Ud). Then it is elementary to verify that = 'Hd(B'^)/'Hd(S'^). It remains 
to show that f{d]a*) > 0. It is easy to see that > Fd+2 > and Fi = l/n < 1 and 
F2 = 1/4 < 1. Thus, 1 — log Fd > 0. Since the digamma function is strictly increasing, 
the expression ip{d/2) — ipi^) ~ log 2 > for d > 4 making f{d; a*) > for d > 4. Direct 
computations show that 

/(I; a*) = [1 + log(7r/8)] /vr, /(2; a*) = 1/4, /(3; a*) = 2 [7 + log(37r/1024)] / (Stt) 

are all positive, which completes the proof of the lower bound. □ 

We need the following auxiliary results for the upper bound in Proposition [3] . Let 
C(x, p) denote the spherical cap {y G S'^ : |y — x| < p} = {y G S'^ : (y, x) > 1 — p^/2}. 

Lemma 5. Let d > 1. Then for x G S'^ and < p < 2, the normalized surface area 
measure of the spherical cap C{'x,p) is given by 

where the Gauss hypergeometric function is a polynomial if d is even and reduces to 1 if 
d = 2. 

Proof. Using the definition of the spherical cap, the Funk-Hecke formula (see pi]), and the 
substitution t = 1 — (p^/2)u, the surface area of a spherical cap can be written in terms of 
a hypergeometric function as follows 

a,(C(x,p))= / dad = ^ r {l-t'f'-'dt 

= Y—P' fu'"-' il-uY-' [1 - (pV4) uY^'-'du 
^ Jo 

_ 1 Ud-^ d nd/2)m A - d/2, d/2_ o . 

-2^P V{l + d/2) l + d/2 

from which the result follows using properties of the Gamma function. □ 
Lemma 6. Let m he a positive integer. For 2; G C \ { — 1, —2, —3, . . . } there holds 

E W^%r [^(^ + ^) - + = 

^-^ k\(m — k]\ mlm 



Proof. Let fm{x) denote the sum for real x > 0. Using the integral representation [36 
Eq. 2.2.4(20)] 

/ dt = V(/3 + l)-^(a + l), Rea,Re/3>-l, 

Jo i- — t 



we obtain 

»1 1 J.X-1 



/m(x) = ^ gm{x,t)dt, g^{x,t):=J2 
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where the function Qmi^, t) can be expressed as regularized Gauss hypergeometric functions 



ml 



1 + X — m' 



1 — m, 1 + X. 
X — m ' 



t). 



Substituting the series expansion of gm{x,t) we have (for x — m not a negative integer) 



r(l + x)!^' (l-m),(l + x), 



-ly 
-ly 

-ly 

-ly 



m\ ^-^ r(l + X — m + /c)/c! Jq 
r(l + x):^' (l-m),(l + x), 1 



- t'^+^-i) dt 



r(l + x):^' (l-m),(l + x). 



m! r(l + X — m + /c)/c! + 1 



■1)' 



E 



m! r(l + X — m + k)k\ k + x 



^ m 

i y 



-m)^T{k + x) 



m\(—m) r(x — m + k)k\ 
1 r(x) ^ 



(1 - m)^ r(A; + x) 
m! r(l + X — m + /;:)/;:! 

m— 1 



■m 1 ^ 



E 



m!(— m) r(x — m) (x — m)j^k\ 

. 1 r(x) 

m!(— m) r(x — m) 



-(-1)' 



r(x) 



(1 -m)^(x)^ 



ml r(x + 1 — m) (x + 1 — ra)^k\ 

( T)m ^ ^(^) p / l-m,x . 1 
^ ^ m!r(x + l-m)' 'U + l-m'^ 



Using the Chu-Vandermonde Identity [li Eq. 15.4.24] 2F1 ^'"'^^ 
obtain 



-ly 

-1)' 
-1)' 



r(x) 



m!(— m) r(x — m) 

. 1 r(x) 

m!(— m) r(x — m) 

. 1 r(x) 



-m] 



_(x - m)^ 



-ly 



1 (l-x)^l 



n,b.,^ = (c-6)„/(c)„, we 

r(x) (1-^)^-1 

m! r(x + 1 — m) (x + 1 — m)^_^ 
1 



1 



■m\mT{x — m) m\m{x)_^ 

where in the last line we used properties of the Pochhammer symbol (see, for example, jSTf 
Appendix II. 2]. Since fm{x) is, in fact, analytic in C with poles at negative integers due 
to the digamma function (the singularity at can be removed), the identity 

(1 - ^In. 



ml m 



can be extended to C \ { — 1, —2, —3, . . . } by analytic continuation. 
Lemma 7. Let d > 1. For x G S"' and < p < 2 



□ 



|x-y| '^dcrd(y) 



§<*\C(x,p) 



-\ogp)-l^ [^{d/2) - ^(1) - 2 log 2] 



1 (^d- 



-E 



2 ujd ^ — ' m\ m 

m=l 



2 uJd 

(1 - d/2)^ (p^ ^rn 

2, 
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where the series terminates after finitely many terms if d is even. 

Proof. Writing |x — y l^'' = [2(1 — t)Y'^^'^ with t = (x, y) and using the substitution 1 + t = 
2(1 — p^/4)m the integral can be expressed as a hypergeometric function as in the proof of 
Lemma [U 

/ |x - yr'^da,(y) = ^ [2 (1 - t)]-'^' (l - t^f'-'dt 

Js''\C(x.p) J-1 



2''/' j''"''^ (1 _ (1 + tY'^-' d t 

(1 - p'/4f f u"^-' (1 - uf-' [1 - (1 - pV4) u] ' d 
Jo 

i (1 - p^Af Tid/2) r(l) 1 - P'/^) ■ 



2 ujd 



u 



2 uJd 

The hnear transformation [Ij Eq. 15.8.10] apphed to the hypergeometric function above 
gives 

|x-y|~'^dord(y) 

n'i\C{yL,p) 

4^ (1 - pV^f' f: ^^^^(^ {Pm'' [2 log(p/2) + i.{k + d/2) - i^ik + 1)] 

k=0 



(1 - ,74)- { 2 (.0, E ^ (0^' + E ^ ^ - . 1)1 (f 



fc=0 A;=0 

^'^{-lo,f.) + '^log2-l'^{l~ f.V4f' ± mk + <i/2) - + 1)] (^ 

ojd iOd I oJd r-~. fc! 

fc=0 

The binomial expansion of the factor (1 — is absolutely convergent for < p < 2 

as is the infinite series above. This gives 

.2, 



2k 



(1 - pV4)^/^ ^ [^(^ + ^/2) - + 1)] ( 

k=0 

- E E ^ + - ii>^ + 1)1 ^ iir^ - E ^..i'') (? 

k=0 n=0 m=0 

where 

fc=0 

Then 6o('^) = '4^{d/2) — and, by Lemma [6| 

&™(rf) = ii^:^, m > 1. 
m! m 



2m 
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Thus, we obtain 

/ |x-yr''do-d(y) = ^:^^(-logp) + ^:^log2 



\C(x,p) 

771=1 



2 Ud 



□ 



Proof of upper hound in Proposition A closer inspection of the second part of the proof 
of Theorem 3 in [27j gives almost (up to a log log N factor) the correct order of the second 
term. Let X*^ = {x^,...,x^} be an optimal d-energy configuration on S*^. For r > 
consider 

N 

D,{ry.=S'\Ci^l,rN-'/''), k = l,...,N, D(r):= fj D,(r). 

k=l 

Since X"^ is a minimal d-energy configuration, for each j = 1, . . . , N the function 
f^.(x):= Yl ^r^' X G S'^ \ (X;^ \ {x*}). 



attains its minimum at x*. By Lemma [7| and ioi p = r N ^/'^ (0 < p < 2) we get 
^,(D(r))f/,(x,)< / f/,(x)da,(x)< I ^ ,, da,(x) 



= (iV - 1) j — {-\og{rN-'")) - [V.(d/2) - ^(1) - 2 log 2] 

\ ^d ^ '^d 

2 ^ 
Hence, 

^,(§^;iV) = ff/.K*) < ^^^^^1^— logiV+ ^ (-logr) - [V.(rf/2) - ^(1) - 21og2] 



1 c.,_i >^ M/2)^ / 1 2\ /riV-V^^^"^ 



2 cjrf (m — 1)! \m (i, 

m=l ^ ' ^ ' 

Subtracting off the dominant term of the asymptotic expansion, we obtain 

f,(S^i N) - l^iV= log N < i ^ i-M) -l/N ^, ^ ^ 1^^^^ (- log r^) 

d ujd d Ud crd{D{r}} d Ud crd{D{r}} 

l^^-lr/^^/o^ 01 01^(^-1) l^d-iNiN-l)^{-d/2)^ fl 2\ /riV-VA^™ 
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Using Lemma [s] in the trivial bound (recall p = rN ^Z'^) 
(36) a,{D{r)) > 1 - Na,{C{^l,rN-'/')) = 1 - ^^r"^ 2F1 " ^Z^, ^/2. 
gives 

[1 - a,{D{r))] logiV < ^'logiV ' ^/^^f 2. 

Choosing r'^ = 1/ log A^, we arrive at the result. □ 

Proof of Proposition \^ We follow the Proof of Proposition [sj now for the kernel Ks{s; t) := 
(2 — 2t + which is positive definite in the sense of Schoenberg |10] with the expansion 

Ke^s-.t) = Yl'^=o^n{s; e) P^{t). (The positivity of the coefficients a„(s;£) can be seen by 
applying Rodrigues formula (see [31]) and integration by parts n times.) We have 

EsiX^) > ao{s; e) - K,{s; 1) N = ao{s; e) - e'^/^ 

where the coefficient ao(s; e) can be expressed in terms of a regularized Gauss hypergeomtric 
function 



ao{s;e) = 



uJd ' Jo \ l + s/A J 

= 2'^— 1^ r(ci/2) r(ci/2) (1 + e/4r^' 2F1 (^M^; . 
ujd \ a 1 + e/A J 

For {s — d)/2 not an integer we can do asymptotic analysis by applying the linear trans- 
formation [H Eq. 15.8.5], that is 

F (s/2.d/2. 1 A [ + p f s/2,1 - d + s/2_ _e_ 

' 'V d 'l + e/4j sm[n{d-s)/2]\r{d-s/2)r{d/2)^ \ l + is~d)/2 ' 4 

[ie/A)/il+e/A)f~'^/\l+e/Ar-^/' fd-s/2,l-s/2 _s 
T{s/2)T{d/2) ' 'V ^-{s-d)/2 ' 4 

which, after simplifications and application of the last transformation in [U Eq.s 15.8.1], 
yields 



2 



r((d-s)/2)sin[7r(d-s)/2] ' 'V ^ + {s-d)/2 
,_i„,u;,_ir(rf/2) TT 



r(s/2)sin[7r(d-s)/2] V4 



^ / l-d/2,d/2__e_\ 
^\l-is-d)/2^ a) 
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Changing to Gauss hypergeometric functions shows that 



ao[s;e) 



1 



ood r(s/2) sin[7r(s - d)/2] T{1 - (s - d)/2) V4 

1 1 71 

- s)/2) r(l - (d - s)/2) sin[7r(rf - s)/2] ^ 



( l-d/2,d/2__e_ 
''\l-{s-d)/T 4 

(s/2,l-d + sl2__e_\ 
^\ l + {s-d)/2 ' 4/ 



Apphcation of the reflection formula for the gamma function p,,, Eq.s 5.5.3] and the sub- 
stitution e/A = a?/('>--'')N~'^/'^ gives 



(37) 

ao(s;£) = 2 



,_,_^Ud-iT{d/2)T{{s-d)/2) 



^2/(d-s)^-2/d 



( l-d/2,d/2_ 
Tis/2) -^i^l-(s-d)/2' 

Note that the first hypergeometric function above is a polynomial if d is even and reduces 
to 1 if c? = 2. Hence, using the series expansion of a hypergeometric function, 



Ssi^'^; N) > f{s, d- a) N^^'/'^ + Vs{^'^) + 0{N^- 



-s/d-2/d 



) as — )■ oo. 



where the function f{s,d;a):=ca — 2 '^a'^^^^ '^^ with (cf. (18)) 



. _ . _ ,d-i-s ^<i-iT{d/2)^{{s~d)/2) _ ^,„i_,r((rf + l)/2) r((s - rf)/2) 



r(s/2) \/^r(s/2) 

has a unique maximum at a* = [2^c{s — d) / s]^^"'^^^'^ with value 



As,d:=f{s, d; a*) = c [2'c {s - d) /s]'^''-^ - 2'' [2'c {s - d) /s]^/^ 



d 



s — d 



,s — d 



s/d 



d 



s — d 



ir((d + i)/2)r(i + (s-rf)/2) 

2 v^r(l + s/2) 



s/d 



□ 



For the proof of Proposition [5] we need the following auxiliary result. 

Lemma 8. Let d > 1 and s > d and {s — d)/2 not an integer. For x G §^ and < p < 2 
we have 



S'i\C{x,p) 



|x-yrMa,(y) = K(S'^) + ^ — 

s — d ood 



'(o/2)'-' F f^-d/2,id-s)/2_p'\ 
~^PI^> l-{s-d)/2 'TJ- 
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Proof. Similar as in the proof of Lemma [7] we obtain 
/ |x - y|-Ma,(y) = ^ [2 (1 - t)]"^/^ (l - t^f'-'dt 

J -I 

^d Jo 

= ^'-'-'^ (1 - pV^f' nd/2) r(i) (1 + 1 - pV4) , 

for {s — d)/2 not a positive integer we can apply the linear transformation [H Eq. 15.8.4], 



,d-i-s^d-i 2 /Ad/2 TT / r(rf/2)r(i) ~ A + (s-d)/2,rf/2 p; 

sm[TT{d-s)/2]\T{d-s/2)T{l)^ \ l + {s-d)/2 '4 



r(rf/2)r(i) 



T{l + {s-d)/2)T{d/2) 



( 2fMd~s)/2 p / rf-s/2,1 _P^\ 



where the first regularized hypergeometric function can be evaluated using [1, Eq. 15.4.6] 
and the second can be transformed by the last linear transformation in [U Eq. 15.8.1], 

^ ^d-i-s^ _ 2/4^/2 ^ f nd/2) . _ 2/^yd/2 

Ud ^ ' ' sin[7r(d-s)/2] 1 r(d-s/2)r(l + (s-d)/2) ^ ' ' 

1 



r(i + (s-d)/2) 

r(rf/2) TT 1 



Wrf r(rf - s/2) sin[7r(rf - s) /2] r(l - (d - s) /2) 

TT 1 rn./o^'^-^ P (l-d/2,{d-s)l2_P^\ 

ujd M<d-s)/2]T{l + {s-d)/2)^P''> l-(s-rf)/2 '4| 

after changing to a non-regularized hypergeometric function and using the reflection for- 
mula for the gamma function, 

_^d-i-s^d-i T{d/2)T{{d-s)/2) , 2'^-a;,_i A - rf/2, (rf - s)/2 

^ r(rf-s/2) ^—d^^P'^^ l-{s-d)/2 'tJ- 



The substitution Ud-i/ujd = T{{d + l)/2)/[y^T{d/2)] (cf. (18)) shows that the first term 
is Vsi^'^) by Eq. The resuh follows. □ 
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Proof of Proposition \^ We follow the proof of the upper bound in Proposition [3] Since 
X'^ is a minimal s-energy configuration, for each j = 1, . . . , the function 



XI — X 



xgS'^\(x;^\{x*}) 



attains its minimum at x*. By Lemma |8j with p = r N-^/^ (0 < p < 2) we get 
a,(D(r))f/,(x,)< / t/,(x)da,(x)< ^ / , . ^ da,(x) 

JDir) fc^fc^^.^D.M |Xfe-X| 



1 ^d- 



Hence 

a,(D(r))^,(§'^;iV)<iV2<'K(S'^ 



1 ^'^-IZ-^AT 



(s - rf)/2 ' 4 



A - d/2, {d - s)/2. r2 /A 
l-(.-rf)/2 'T^ J 



1 uJd- 



S — d LOd 



By relations (19) and (36), we have 



Note, that for = 2 the hypergeometric function in (36) reduces to one and the (9(.)-term 
above disappears. Hence 

;r,d—S 



1 ^d-\ 



s-d Ud 1 - ^ 



The function /i(r) = r*^ — cr'^) (where c = (l/(i)(co'rf_i/u;rf)) has a single minimum in 
the interval (0, 00) at r* = c-^/'^(l - rf/s)^/'^ with value h{r*) = (s/rf)c-^+"/'^(l - d/sf-'/'^, 
where 1 — cr"^ = d/ s > 0. The result follows. □ 



7. Motivations for conjectures 
Motivation for Conjecture^ Suppose Conjecture [3] holds. Proceeding formally, we obtain 



d 
d7 



|Adi 



where 



A(§^; A^) = lim - 

s-s>0+ S 



UJd 



s/d 
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Assuming the limit exists and that A(§^; A^) = o{N) as — t- oo, we have 
(38) 
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^iog(S';Ar) 



ds 



+ oiN), 



M„g(S'')iV2 + CA,(0)^iVlogiV + 



V,,^(S'')N''-2N\ogN + 



d 



^d 

N + o{N), 



N + o{N), 



ViagiE,'')N' --N\ogN + Ciog,dN + o{N) as N oo, 



where (cf. ^ and (10)) we used 



ds 



K(§' 



d 
d^ 



^,^,_, T{id + l)/2)T{{d-s)/2) 
y/^T{d-s/2) 



\og^+lmd)-Hdm] = v,o,{s''), 



and also used the fact that Ca(0) = ^1 holds for any lattice A (cf. |i6]). 

Using Proposition 6] in the appendix and IA2I = (for the hexagonal lattice with 

unit length edges), we obtain 
(39) 

aog.2 = C'A.(0)-^log^ = log(27r)-^-31ogr(l/3)-^log^ = 21og2+^log^+31og^^, 

□ 

Motivation for Conjecture\E^ We first remark that ^^(S'^) has a simple pole at s = d with 



s — d 



+ Ad + 0{\s - d\) as s d, 



where 



and 



a_i,rf:=Res,=dK(§') = -d 



Ad.= lim 



s^d 



s — d 



1 



2 iOd 



i(7-21og2 + V^(rf/2)). 



In addition to Conjecture |3| we further assume that Cs^ behaves, near s = d, as follows 

Cs,d 



where 



6_i,rf/(s -d) + Bd + 0{\s - d\) as s ^d, 



6_i w := -a_i w and i?d:=lim 



[Hdm 



d]]''/'^ s - d 
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Proceeding similarly as before and taking s — )■ d, we have 



[Ham 



s — d 



+ 



C. 



s,d 



s — d 



l'U^(Sd)Y/d s - d 
^^dN' + ^ + ^dN' + A,{S'; AT), 

where we assume the limit Ad{S'^; N) exists and that AdiS"^; N) = o(A^^) as N ^ oo. 

Multiplying both sides of Q with [Hdi^'^)]'^'^ and expanding '^'^ about s = d, 

we obtain that 

(40) 
where 



C. 



s — d 



+ Bd + 0{\s-d\) as s ^ d, 



5rf:=lim 

s^d 



u^d-i 

s — d 



Furthermore, the following relation holds between coefficient of the A^^-term and Bd'. 

1 UJd-l , _ , _ , . . ,„x X 1 Wri-l , Bd 



Cd4.=Ad + Bd = -TT— (7 - 21og2 + ^{d/2)) 
I bJd 



1 Wd-l , , Dd 

-, logWd H , 

d ujd ujd 



In the case d = 2 (with the help of Mathematica), we obtain 



a-1,2 



Res 



s=2 



2-s 



-1/2, 6_i^2 = Res^=2 



s,2 



(An 



«/2 



1/2, 



and 



= lim 



)1-S 



Bo = lim 



s-5>2 



a~i,2 

2-s s-2 

Cs,2 &-l,2 



log 2 



(47r)^/' s - 2 



4 L 



7-log(8v^7r) +^[7i(2/3)-7i(l/3)] 



An 



where 7 is the Euler-Mascheroni constant and 7n(a) is the generalized Stieltjes constant 
appearing as the coefficient of (1 — s)" in the expansion of ({s, a) about s = 1. 

□ 



Proposition 6. 

L_3(0) = 1/3, 
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Appendix A. Auxiliary results 
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Ca,(0) = -1, G,(0)=log 



Proof. Note that one has the following identities (cf. [HI P- 264]) 

d 



(41) 



C(0,a) = (l/2)-a, 



ds 



C(s,a) 



logr(a) - logV27r 



By (|29|) and the last two relations 
,l/3)-C(0,2/3) = ( 

{3- [C(s,l/3)-C(s,2/3)]} 



L_3(0) = C(0, 1/3) - C(0, 2/3) = (1/2) - (1/3) - (1/2) + (2/3) = 1/3 
d^ 



LU(0) 



- log3 L_3(0) + logr(l/3) - logr(2/3). 



S-5>0 



For the second relation for L'_3(0) we used r(2/3) = 27r/[v^r(l/3)]. 
By ^ and C(0) = -1/2, we get 

Ca,(0) = 6 C(0) L_3(0) = 6 (-1/2) (1/3) = -1 



ds 



{6C(s/2)L_3(3/2)} 



= 3C'(0)L_3(0) + 3C(0) L'_3(0). 

) 

C(s, 1), we derive from (41) the special values ({0) = —1/2 and C'(0) 



Since ({s) 
— log \/2n. This completes the proof. 



□ 
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